Measuring the local dark matter density 
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Abstract. We examine systematic problems in determining the local matter density from the 
vertical motion of stars, i.e. the "Oort limit". Using collisionless simulations and a Monte Carlo 
Markov Chain technique, we determine the data quality required to detect local dark matter at its 
expected density. We find that systematic errors are more important than observational errors and 
apply our technique to Hipparcos data to reassign realistic error bars to the local dark matter density. 
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Introduction. Determining the local density of matter is an old problem [2] with 
a renewed importance for direct dark matter (DM) searches. The local mass density 
determined from measurements of stellar motions is 0.102±0.01M Q pc _3 [3, 4], limited 
by the sample size of stars with good distance estimates and proper motions. The quoted 
uncertainties are of the same order as the DM density in the standard dark halo model 
for the Galaxy. Upcoming surveys promise a dramatic improvement in the number and 
precision of astrometric, photometric and spectroscopic measurements. 
With this in mind, we revisit the Oort limit [1]. Following the work of [3, 5, 6], we 
determine the local DM density by comparing the predicted density v(z) found from the 
vertical velocity distribution at the plane /(wq) of a stellar tracer, (equation 1), 



with the measured stellar density fall-off in cylindrical patch at the Solar Neighborhood. 
As in [3, 5], we calculate the potential by modeling the local visible mass as a super- 
position of isothermal components, and the local dark mass as a constant density term. 
We then integrate the Poisson and Boltzmann equations simultaneously, neglecting the 
cross terms of the latter. These approximations are valid as long as the potential can be 
separated into its vertical and radial component, <t>(R, z) = 4>(z) + <!>(/?) , the local visible 
matter is made up of nearly isothermal species and the disk scale height is much smaller 
than the disk and the halo scale lengths: Zd <C r^ <C R^. 

Analysis and results. We apply this method to a high resolution (30 million disc 
particles) collisionless simulation of a MW-like galaxy, whose initial condition were 
produced using the galaxy models by Widrow and Dubinski [7], and fit the density fall- 
off with a Monte Carlo Markov Chain technique to find the best value of of the local 
stellar and DM densities (p s , and Pdm)- First, we take an early snapshot of the simulation 
(~50Myrs) when the disc is still axisymmetric to produce a "super-patch" containing a 
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FIGURE 1. Recovered star (upper rectangles) and DM (lower rectangles) densities for the axisymmetric 
(left) and inhomogeneous (right) discs in four volumes around the disc. The shaded regions represent the 
density of the MCMC models and the crosses the best % 2 values. 

large number of particles at the solar radius, comparable to what we will get with the 
next generation surveys (e.g. GAIA). Then we use the evolved simulation (~0.7Gyrs), 
which has a bar and spiral waves, to evaluate the impact of disc inhomogeneities. 
Taking stars in larger volumes around the Sun improves the statistics, but it's important 
to evaluate how far we can go in z and R before breaking the method's assumptions. 
From the analysis of the axisymmetric disc, using boxes with different radial and vertical 
sizes, we can recover the local p s and Pdm provided we use a local volume to compute 
the velocity distribution with \z\ < 75pc, R < 300pc and we fit the density fall-off up to 
3-4 times the disc scale height Zd (i.e. z > 0.75kpc). This latter is required to break a 
strong degeneracy between p s and Pom due to the potential similarities at low z ■ 
In the axisymmetric disc, we recover the local density correctly within la (of the 
MCMC model distribution) as long as we have a good sampling (e.g. GAIA data) and the 
systematic errors are of the order of the observational uncertainties and of the expected 
amount of DM in the standard dark halo model for the Galaxy (~ O.OlM pc~ 3 ), as 
shown in figure 1 . 

The situation is different in presence of disc inhomogeneities: in this case the errors are 
larger and the method seems to systematically overestimate the stellar density (figure 
1 right). This may owe to the cross terms of the velocity dispersion tensor in the 
inhomogeneous disc, that are neglected in the method. 

The "super patches" have statistics that are comparable to the future GAIA sample, 
so we want to introduce realistic errors in the velocities and positions of the stars to 
quantify their effects on the recovered density. We find that these errors don't change 
the recovered values for the density, but merely increase the % 2 of the fit. This means 
that systematic errors have a dominant effect. For this reason, it is important to calculate 
the local matter density for the available data, using our MCMC technique, to assign 
realistic errors to the DM density measure. This will be the next step of our analysis. 
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